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Abstract 

In kaon decay, electromagnetic radiative corrections can generate shifts 
in the apparent AI — 3/2 amplitude of order olAqI A 2 ~ 22a. In or- 
der to know the true AI — 3/2 amplitude for comparison with lattice 
calculations and phenomenology, one needs to subtract off this electro- 
magnetic effect. We provide a careful estimate of the leading electro- 
magnetic shift in the chiral expansion of the amplitude, which shows 
that it is smaller than naive expectations, with a fractional shift of 
5A { ° m) /A 2 = -0.016 ±0.01. 



1 Introduction 



The AI = 1/2 enhancement of nonleptonic kaon decays is a well-known 
feature which has still not been completely explained. Most existing lattice 
calculations find a AI = 3/2 amplitude which is a factor of two larger than 
the experimental value.[] This AI = 3/2 amplitude also enters present phe- 
nomenology through the chiral determination of the Bk parameter which 
is part of the Standard Model prediction of CP violation. The chiral cal- 
culation [H] (which relates Bk to the experimental AI = 3/2 amplitude) 
disagrees with the quenched lattice calculation for B^ also by about a fac- 
tor of two. At present we do not know if the disagreement in Bk is due to 
the failure of the lattice approach to reproduce the experimental amplitude 
or if there are large chiral corrections responsible for the difference. 

In both these applications, we need to know the true AI = 3/2 ef- 
fect. Since the AI = 1/2 amplitude is so much larger, it is possible that 
electromagnetic corrections to it may simulate an effect which is similar to 
the small AI = 3/2 amplitude. These electromagnetic corrections enter 
at order aAo = Ao/137, which can be comparable to a sizeable portion of 
the AI = 3/2 amplitude A2 — Aq/22. The possibility then emerges that 
the relevant experimental amplitude (without electromagnetism) could differ 
significantly from that presently being used in phenomenology. 

Although the literature on this subject extends over many years, e.g. see 
Refs. H ||, ||, g|, the issue has not yet received a definitive treatment. It 
is our aim in this paper to provide an analysis using the most up-to-date 
tools which will yield a reliable estimate of this effect. In a longer paper, 
we will present a more comprehensive analysis of electromagnetic radiative 
corrections in kaon decays. The full system, particularly the decay K$ — ► 
7r + 7r~, brings in several additional complications, such as the Coulomb effect 
on the final state, the violations of Watson's theorem from the mixing of final 
states and the induced AI = 5/2 effect. However, the decay K + — > 7r + 7r° 
is particularly simple and can by itself be used to answer the question that 
we have raised above. Among our results in this paper are: 

1. The long distance portion of the leading chiral result satisfies the re- 
lation AihD = — 2g8(3 m n)LD / F% where (dm^)^ is the long distance 
part of the pion electromagnetic mass difference and the other con- 
stants are defined below. 

1 A recent lattice calculation is more promising in this regard. Jl]] 
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Figure 1: Combined weak-electromagnetic transitions of quarks. 

2. This leading long distance effect is canceled exactly by the effect of 
pion mass differences in the usual weak amplitude. 

3. There are, however, residual effects coming from intermediate energies 
and the electroweak penguin operator. Although we allow generous 
uncertainties associated with intermediate energies, it is clear that the 
net residual effect is quite small. 

2 Chiral Analysis 

Chiral symmetry provides the framework for structuring our calculation. 
We first replace the calculation of the decay amplitude A + q(pk,P- k + iP-k ) 
with the simpler K + -to-7r + matrix element by taking the soft pion limit, 

A +0 (p,p,0) = -^-M(p) , (f) 

where 

out (7r+(p / )|^ + (p))in = ^(2^) 4 ^ 4 )(p'-p) M{p) . (2) 

For most of our analysis, we shall work with the leading chiral component 
M(0) 

M{p) = M(0) + 0{p 2 ) + ... , (3) 
and add in (small) 0{p 2 ) contributions at the end. 

2.1 The Chiral and Electromagnetic-penguin Components 

At the quark level, electromagnetic corrections to the weak transition s+u — ► 
d + u fall into two distinct classes (of Fig. 1), which we call the chiral (CH) 
and electromagnetic-penguin (EMP) components, 

M(0) =Mch + -Memp • (4) 
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The amplitude AfcH is associated with the long-range and intermediate- 
range contributions of the process in Fig. 1(a) (along with all other diagrams 
in which a photon is exchanged between the quark legs). We note that the 
short distance part of such transitions leads merely to an overall shift in the 
strength of the nonleptonic interaction. Using the procedure described in 
Ref. this is equivalent to an effective Fermi constant defined as 



where p, is an energy scale lying in the region where perturbation theory is 
valid and is the Fermi constant measured in muon decay. This shift does 
not lead to mixing of isospin amplitudes and is irrelevant for the purposes 
of this paper. Calculation of the long and intermediate range contributions 
to AIch is carried out in Sect. 3.1 and Sect. 3.2. 

The amplitude AIemp corresponding to the electromagnetic penguin op- 
erator of Fig. 1(b) will have both long-distance and short-distance compo- 
nents. These are described in Sect. 3.3. Determination of the full amplitude 
M. is carried out in Sect. 3.4, with special attention paid to the relative 
phase between AIch and AIemp and to the matching of long and short 
distances. 

2.2 Chiral Lagrangians 

We now introduce some useful chiral lagrangians. The weak interactions 
involve left-handed currents only and the nonleptonic hamiltonian has an 
octet and 27-plet component. The lowest-order weak lagrangian for the 
dominant octet portion involves two derivatives, 



with la, | ~ 7.8 • 1(T 8 F%. 

Electromagnetic corrections involve both left-handed and right-handed 
effects, and can lead to lagrangians which do not involve derivatives. For 
example, one of the effects of electromagnetism is to shift the charged pion 
masses, an effect described at lowest order by the lagrangian 



The parameter g ems is fixed from the pion electromagnetic mass splitting, 






(6) 




(7) 



(8) 



7T 
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There is also the lagrangian which describes the leading electromagnetic 
correction to the weak interactions to leading chiral order .0 

£emw =fawTr (\&UQU^) , (9) 

where 

jcmw is an a, priori unknown coupling constant. Knowledge of fl'emw 
is equivalent to that of the matrix element A4(0) as the two are related in 
the chiral limit by 

M(0) = 5emw • (10) 

In Sect. 3, we present a detailed calculation of M(0) (and thus of j emw ) 
and as a consequence reveal an approximate numerical relationship between 
5emw and <? e ms- Then in Sect. 4, we turn to the full A + o amplitude, including 
also the effects that arise at the next chiral order (0(p 2 )) from the photon 
loop calculation. 

3 Calculation of the leading chiral amplitude 

To fully calculate the relevant amplitude, we need to consider contributions 
from all scales. We will recognize three regions of the virtual photon mo- 
mentum: 

1. very low energies Q 2 < A 2 with A ~ m p , 

2. high energies with Q 2 > ii 2 (n ~ 1.5 -> 2.5 GeV), 

3. intermediate energies between these two regions. 

In the lowest energy regime, we may use chiral techniques to obtain the 
leading effect. At high energies, the short distance analysis of QCD will be 
employed. The most important ingredient of the treatment of intermediate 
energies is the requirement of matching these two regions. This will be 
modelled on physics which is reliably known in the case of electromagnetic 
mass shifts. 

2 We also make use of lagrangians which couple pions and kaons to resonances, as in 
Figs. 2(b), 3(b). 
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(a) 




Figure 2: Electromagnetic corrections to the K + -to-ir + transition. 



3.1 Long Distance Component of Aicu 

The very long distance component can be calculated from the combined 
chiral lagrangian of the strong weak and electromagnetic interactions in the 
chiral limit. In the diagrams of Fig. 2(a), one finds after Wick rotation a 
matrix element 

where A represents the upper end of the low-energy region. We note that 
a similar calculation of the electromagnetic mass shift of the charged pion, 
Fig. 3(a), yields 

6 <L = If > ■ < 12 > 

The similarity of the two can be motivated by the fact that in the former 
calculation the weak vertex in the loop introduces a factor of I 2 (I is the 
loop momentum) which compensates one of the two propagators, yielding 
an effect similar to that of Fig. 3(a). At this stage, it is a curiosity to note 
that the choice of A 2 = m? p provides an accurate description of the pion 
mass difference. However, we will see below that this is not an accident - 
that reliably known physics cuts off the integral above the rho mass. For our 
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(a) (b) 
Figure 3: Pion electromagnetic mass shift. 



purposes at this stage, this similarity is the first indication of the relation 

M LD = -2^(5ml) LD (13) 

or (<?emw)LD = ~98^ r n 2 LD- This could be derived somewhat more formally 
by using a rotation to the basis where the kinetic energy matrix is diagonal- 
ized. H The application of long distance electromagnetic mass shifts to the 
rotated basis is equivalent to the above relation in the non-rotated basis. 



3.2 Intermediate Energy Component of .M CH 



A prototype for dealing with the intermediate energies is the pion electro- 
magnetic mass difference. The accumulated wisdom of many studies has 
given us an accurate guide to the physics of this process. A rigorous ap- 
proach would involve the sum rule of Das et al jR]] , in which Sm 2 is expressed 
in terms of the difference of the experimental vector and axialvector spectral 
functions (py— pa)(s). This has been analysed successfully using experimen- 



tal data and QCD constraints. [11| A simplified expression that captures the 
essential physics is obtained upon saturating py and pa respectively with 
the vector p resonance and the axialvector resonance a±. This yields 
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(14) 



The second form is found when the resonance couplings and masses satisfy 
the Weinberg sum rules, which are required in order to obtain the right high 
energy behavior. The long-distance amplitude given in Eq. (|l^) has been 
softened at values of Q 2 above the meson masses so that m p and m ai act as 
the effective cutoff for the integral. This result is equally well reproduced 
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by introducing resonance couplings to the effective lagrangian and imposing 
the Weinberg sum rules on the masses and couplings. This involves the 
diagrams of Fig. 3(b). 

What is the analogous statement for Mcu = -Mld + -Mint + -Msd? 
First, we recall from the discussion surrounding Eq. (|5|) that, rather than 
contributing to the isospin mixing effect, the dominant effect of the short 
distance (SD) component is to renormalize the Fermi constant. The full 
chiral amplitude thus experiences important contributions only from long 
and intermediate distance effects and must vanish in the short distance 
region. Combining results from the previous sections, our form for the long 
and intermediate distance regions is 

B V Q 2 _ B A Q 2 Cm\Q 2 
Q 2 + ra 2 / Q 2 + m\ (Q 2 + m 2 A ) 2 \ ' 

(15) 

where the Q 2 -integral for Aicn is seen to be effectively cut off at some scale 
fj?. The quantities By, B A and C in the above contain couplings from 
the weak interaction resonance lagrangians. It is of course possible to 
model these couplings, and we have done so. However, it is more to the 
point to implement the constraint (discussed earlier) that TWch receive no 
mixing contribution from the high-Q 2 region. Thus these constants must 
be constrained such that the ampliude vanish at the matching to the short 
distance region. For this to occur, we require that the large Q 2 limit of the 
integrand (Q 2 3> ray a) approach zero rather than a constant, i.e. 

B v + B A = 1 . (16) 

We further constrain this amplidude by choosing the matching scale Q = [i 
at which the amplitude vanish. If ray and ra A were equal this constraint 
would determine unknown C in the integrand of Eq. (|l5| ) . In our numerical 
study, we choose \x to lie between 1.5 GeV and 2.5 GeV and treat the 
resulting variation as one of the uncertainties of the calculation. Of course, 
the difference between ray and tha leads to a slight further uncertainty. To 
model this effect, we have explored models for the resonance couplings 
— which weight the axialvector and vector resonances differently. It is found 
that this uncertainty is smaller than that associated with variation of the 
matching scale 



3ag 8 2 



7 







u, c f 


) % c 


/ 


> 7 


u 


u 



Figure 4: Electromagnetic penguin and the i^ + -to-7r + transition. 
3.3 Determination of .Memp 

As we shall show in the following, the electromagnetic penguin (EMP) op- 
erator of Fig. 1(b) gives rise to contributions over all distance scales. It is, 
however, the sole source of meaningful short distance effects in our calcu- 
lation of electromagnetic corrections to A2 ■ If the numerically tiny i-quark 
contribution is omitted, the EMP hamiltonian takes the form, 



-iHemp = G 



d 4 g J^(g,/2) 
(2vr) 4 q 2 + ie 



x / d A y T[s(0)M l +^)d(0)q(y)QlMy)\ 



(17) 



where q = u, d, s is a light-quark field, G = 2e 2 G F V* s V ud /(3V2), and 
I IJ,v {q^fi) represents the effect of the quark-antiquark loop in the EMP op- 
erator. In evaluating I^ u (q,fl), it is understood that at the lower end, the 



loop momentum is cut off at scale fi (c/ see Eq. (21)). The dependence on 
ft is logarithmic and thus quite weak. 



It can be shown [14| that in the chiral limit the K + -to-ir + matrix element 
of "Hemp ( C /Fig. 4) is expressible as 



lim (7T + (p)\T [s(0) 7/1 (l + K)d(0)q(y)Qj„q(y)] \K+(p)) 
p=0 



? r d 4 k / \ 

'■p2 J (2^)4 elh ' V ~ fe V) Pvs - n A3 ] (k 2 



-(18) 



where IIv3, IIa3 are the isospin vector and axialvector correlators. In the 
chiral limit, the K + -to-n + matrix element of the EMP operator thus be- 
comes 



lim (ir + (p)\H E Mp\K + (p)) 
p=o 



2G 

f2 



A 7^(g,/2) 
(2vr) 4 q 2 +ie 



[n V3 - n A3 ] te 2 ) -(19) 
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Figure 5: Schematic of matching. 



This expression describes the EMP effect over all scales of the virtual photon 
(euclidean) momentum, 

3GMyM\ Q 2 



where we have expressed the Q 2 dependence of the correlators in terms of 
vector and axialvector pole terms, an approximation we know to be valid to 
within a few per cent. An explicit form for the EMP integral I(Q 2 ,£i 2 ) is 



I(Q 2 ,J1 2 ) = f 1 dx x(l 
Jo 



X 



\i + m 2 + Q x(l — x) 
p? + m 2 + Q 2 x(l — x) 



m 2 + Q 2 x(l-x) m 2 l + Q 2 x(l-x) 



/j, 2 + m 2 + Q 2 x(l — x) /j, 2 + m 2 + Q 2 x(l 



(21) 



where m c and m u are the c-quark and n-quark masses. The latter vanishes 
in the chiral limit. 



3.4 Matching 

The final step in determining A4(0) is to add together the components Men 
and .Memp- This is displayed schematically in Fig. 5 which depicts the inte- 
grands in the Q 2 integrals and indicates that the short-distance contribution 
is numerically much smaller than the long-distance contribution. 



It would appear that the calculation must contain an ambiguity arising 
from ignorance of the relative phase between Mch and A^emp or equivalen- 
tally of the sign of $8- However, we can infer that g$ < from the following 
argument. From the AI =1/2 chiral lagrangian of Eq. (|6|), we have for the 
K° — > (27r)/ = o amplitude, 

A = (m 2 K - ml) g 8 . (22) 

Alternatively, we can obtain Aq using the effective lagrangian of the Stan- 
dard Model 

C AS=1 = -^V ud V: s J2 ctMOi , (23) 
together with the vacuum saturation approximation (VSA), 

((2ir) I=0 \Oi\K°) = B^((2tt) i=0 \O 1 \K ) ysa , (24) 



to write 



with 



A = -iG F F w V ud V: s (m 2 K - ml) gf , (25) 



.eff 



- Ci b[ 0) + -C 2 Bf ] + -C3S3 + C4B4 

d y ~6 



(26) 



In the above, L5 is a coefficient in the Gasser-Leutwyler 0(p ) chiral la- 
grangian, the superscripts on b[ , B^ signify 1 = for the final state 2tt 
pair, and we refer the reader to Ref. [15] for further details. Since (with the 
exception of B3) the {Bi} are all positive, the {\ci\} have been determined 
at NLO and specifically ci,cq < 0, we conclude with reasonable certainty 
that g$ < 0. 

The analysis described throughout this section then leads to the following 
value, expressed in units of g%5ml, for the coupling g cwp of Eq. (||), 

-0.62 ±0.19 . (27) 



g%bm\ 



The uncertainty arises almost entirely from variation of the parameter \x (we 
have set p, = 1.5 GeV in this determination). 
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4 The K 



7r + 7r transition 



We now turn to the physical K + — > 7r + 7r° transition. This receives con- 



tributions from the true AI = 3/2 interaction {A 2 



(true) > 



electromagnetic 



corrections {5A^ ) and isospin-breaking effects (SA^ b ) 



, phys 



— „ 



.4 



(true) 



+ ,54 cm) + M 



(iso— brk) 



(28) 



where 



5A 



(em) 



1 58 



5m£ - <5 A 



(h-o) 



(29) 



The first term in Eq. fl29[ ) was the subject of the analysis in Sect. Ill and 
has been discussed in great detail. The next term has a somewhat subtle 
origin. The contribution from the AI =1/2 weak interaction £$ of Eq. (||) 
to the K + — ► 7r + 7T° amplitude is proportional to g%{jp\+ ~ P^o) = g&Sm^. 
This is ordinarily discarded in calculations in which isospin is conserved and 
Sm 2 = 0. It cannot, however, be neglected in the present context. Perhaps 
the most interesting feature of our result is an approximate cancellation 



between the first and second terms of Eq. (29). To make this explicit we 
recall Eq. (|i~3[) to write 



-g$8m„ + Sg 



(30) 



where Sg emw arises from the sum of the intermediate-range part of the chiral 
contribution and the EMP contribution. Our estimate implies 



$g c 



g 8 5ml 



0.38 ±0.19 



(31) 



for this quantity. 

Finally, the contribution SA^ °^ in Eq. (|29|) represents electromagnetic 
corrections of higher order in the chiral expansion which vanish in the chiral 
limit. The Feynman diagrams for the photonic corrections to the K + — > 
7r + 7T° amplitude also generate effects at order (e 2 p 2 ), and we find 



5 A. 



(h-o) _ 3«g 8 
" 4vrF3 



mL In 



A 2 



rat 



° 2 

+ 2 m ^ + 



(32) 



In 5 A 



(h-o) 



there is a residual dependence on the cutoff A. However, since 



it enters only logarithmically and is multiplied by a small factor of m\ it is 
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inconsequential for the final answer. We have simply set A 2 = m 2 in this 
contribution. In addition, at this order in the chiral expansion one also needs 
to include meson loop effects, involving both the effects of 5m 2 and the loops 
proportional to g e mw This can be done in chiral perturbation theory. Our 
evaluation of g emw becomes an input in that calculation and the photon loop 
result of Eq. (^) is relevant for the determination of the chiral constants at 
order (e 2 p 2 ). At this stage, we will include only the effects of Eq. (|32"|) and 



reserve a full calculation at next order for a future publication [17]. 

Overall, the net result of our calculation is a shift due to electromagnetic 
effects in the apparent A2 amplitude ranging over 0.6 — ► 2.6% depending on 
how the matching is carried out. Taking the mean value, we obtain 



5A 



(cm) 
2 



A-2 



-0.016 ±0.01 , (33) 



We do not calculate the contribution 8A% s ° brk ^ in Eq. ( p8|) which arises 
from the mixing between tt° and 77,?/, 

^(iso-brk) = §A ^- V + §A n°-rf ^ 

This effect is primarily due to quark mass differences and has already been 
analyzed in Ref. [||. The result cited there isf] 

— 2- = — \ + — \ -0.14 + 0.21 = 0.35. 35 

Ao Ao An 



5 Conclusions 

In the K — ► 2tt amplitudes, the ratio of AI = 1/2 and AI = 3/2 amplitudes 
is about 22. This suggests that electromagnetic corrections to the former can 
lead to contributions to the latter of order 22/137 or around 20%. Indeed, we 
have found individual contributions of order 10% to occur. If added up, these 
electromagnetic corrections to the K + — > 7r + 7r° amplitude would contribute 
at the 20% level. However, there turns out to be a significant cancellation 
which greatly weakens the effect. The realization of this cancellation requires 
a consistent application of electromagnetic effects to both the pion masses 
and the weak amplitude. Including all scales in the electromagnetic shifts 
leads to this cancellation only being partial. Due to this cancellation, we 

3 Note that Eq. (IX-3.21) of Ref. ^ contains a minus sign typo. 
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have assigned a generous fractional uncertainty to the final answer. However 
in absolute terms, the overall electromagnetic effect that we have calculated 
at this order is only a small part of the experimental value of Ai- Given the 
knowledge of the long and short distance components of the amplitude, our 
confidence in this result is quite strong. 
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